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Analytic framework and ambipolar transients 
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Up to now the basic theoretical description of charge extraction by linearly increasing voltage (CELIV) is 
solved for a low conductivity approximation only. Here we present the full analytical solution, thus generalize 
the theoretical framework for this method. We compare the analytical solution and the approximated theory, 
showing that especially for typical organic solar cell materials the latter approach has a very limited validity 
Photo-CELIV measurements on poly(3-hexyl thiophene-2,5-diyl):[6,6]-phenyl-C6i butyric acid methyl ester 
based solar cells were then evaluated by fitting the current transients to the analytical solution. We found that 
the fit results are in a very good agreement with the experimental observations, if ambipolar transport is taken 
into account, the origin of which we will discuss. Furthermore we present parametric equations for the mobility 
and the charge carrier density, which can be applied over the entire experimental range of parameters. 



I. INTRODUCTION 

Bulk heteroj unction solar cells use a phase separated 
blend of an electron accepting and electron donating ma- 
terial — e.g., poly(3-hexyl thiophene-2,5-diyl) (P3HT) and 
[6,6]-phenyl-C6i butyric acid methyl ester PCBM — as active 
layer [1]. In this nano-scale blended film the photo-generated 
charges are separated at the donor-acceptor interface and col- 
lected at the electrodes under short circuit conditions. State- 
of-the-art polymer-based solar cells provide high yields for 
collected charges with respect to the incident photons [2] and 
have reached a power conversion efficiency (P.C.E.) up to 
7.9% under AM1.5 (lOOmW/cm^) illumination [3, 4]. 

The knowledge about the physical processes, such as re- 
combination and charge carrier transport, and their impact on 
the charge collection in organic solar cells is crucially impor- 
tant for an optimization of the P.C.E and has therefore been in- 
tensively debated in literature [5, 6]. The recombination and 
the charge transport are governed by the charge carrier life- 
time X and the charge carrier mobility ^, respectively. Various 
techniques have been used to study charge carrier dynamics 
in these systems, e.g. time-of-flight photo-conductivity [7] 
for the transport or transient absorption [8-10] for the recom- 
bination. Another approach to measure charge carrier mo- 
bility and recombination was introduced about ten years ago 
by G. Juska etal. [11], called CELIV, charge carrier extrac- 
tion by linearly increasing voltage. Due to its ability to mea- 
sure these two parameters simultaneously this technique has 
attracted much interest in the organic semiconductor research. 
However, the theory for calculating the current response due 
to the linearly increasing voltage has only been presented for 
the simplified cases of low and high conductivity regimes, re- 
spectively [1 1]. As we will see later, to our knowledge neither 
the low conductivity approximation = ££o / enjd ^ ttr (Xo is 
the relaxation time, ttr is the time at which all free charge car- 
riers are extracted) nor the high conductivity approximation 
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<C ttr is valid for state-of-the-art organic solar cell materi- 
als, such as P3HT and PCBM. 

In this paper we present the full analytical framework 
for the CELIV method which allows us to computationally 
evaluate the mobility fi and the charge carrier density n of 
experimental measurements. In general, a closed analytical 
expression for the relation ju{tmax) cannot be derived, as for 
the high or the low conductivity approximation. Therefore 
we propose a parametric equation for the mobility eval- 
uation from CELIV experiments which is valid over the 
entire experimental range of conductivities. An equivalent 
expression is derived for the charge carrier density n{tynax)- 
The mobility equation is tested against parametric mobility 
equations known from the literature [5, 11-13]. Furthermore 
we evaluate experimental photo-current transients by fitting 
the data directly to the CELIV framework. 



II. CELIV THEORY 
A. Method summary 

A schematic illustration of the CELIV technique is given in 
Fig. 1. In this method a linearly increasing voltage V{t) =A't, 
where A' is the slope of the applied voltage pulse, is used to 
extract equilibrium charge carriers with density n and mobil- 
ity /i from a film with a certain dielectric permittivity 8 and 
thickness d. The whole device is represented as a capacitor 
with the film between two electrodes at x = (blocking con- 
tact) and x = d. We here summarize the important parts of the 
CELIV theory [11] for comprehension and later discussions. 

It is assumed that one charge carrier is much more mobile 
than the other (unipolar transport), that the electrode dimen- 
sions are much larger than the device thickness and the charge 
carrier density of free charge carriers is n. This yields the fol- 
lowing charge carrier distribution p(x,r) at time t 
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FIG. 1 . (colour online) Schematic illustration of the CELIV method, 
(a) Process of charge extraction and the band diagram in the active 
material sandwiched between two electrodes at a certain time. The 
region < x < l{t) where all electrons have already been extracted 
is depicted by the shaded area, (b) Scheme of the voltage input and 
(c) the current output. The voltage pulse with slope A' is applied 
in reverse bias. This yields a characteristic current density response 
with a capacitive offset Jq and current density due to the drift of the 
free charge carriers with a maximum current of A j. 



This expression is related to the time dependent extraction 
depth l{t), where < / (^) < J. At a certain time t in the region 
< X < l{t) all electrons have been extracted and the region 
is charged positively. 

Applying the increasing voltage V{t) with the condition 
RC <C t the total current density j{t) in the external circuit 
due to the redistribution of the charge carriers (electric field) 
is 
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R and C are the device resistance and the device capaci- 
tance, respectively. 70 is the differentiating initial step of the 
RC-circuit, p(^) is the density of free charges in the device 
and becomes zero when l{t) = d. The last term in brackets 
in Eqn. (2a) describes the drift of the free charge carriers due 
to the external field and due to the electric field caused by the 
charge distribution in the sample. In the simple case the drift 
due to the latter electric field is neglected. For the calcula- 
tion of the current transient j{t) the extraction depth /(^) is the 
crucial parameter. l{t) can be expressed as a Ricatti-type first 



order differential equation 
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and has up to now either been solved only numerically [13] or 
analytically for a low (Xq ^ ttr) and a high (Xq ^ ttr) conduc- 
tivity approximation [11, 14], respectively. Note that in the 
low conductivity approach Eqn. (3) simplifies to dl{t)/dt = 
juA't/d and thus the extraction depth l{t) becomes 
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From Eqn. (6) the transit time ttr — corresponding tol{t) =d 
when all free charge carriers are extracted — can be defined 
as 
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In the next section we introduce the analytical solution for the 
Ricatti-type definition of the extraction depth l{t) (Eqn. (3)). 
By doing this, the CELIV theory is generalized and can be 
compared with the low conductivity approach. This way we 
can prove that the low conductivity approximation has a very 
limited validity for state-of-the-art disordered organic semi- 
conductors used in organic electronic devices such as organic 
solar cells. 



B. Solving the Ricatti equation 



For solving Eqn. (3) directly we need to use the substitution 
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l{t)-L{t). 
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Eqn. (8) is again a differential equation for the substituent 
L{t). Inserting Eqn. (8) in Eqn. (3) transforms the latter equa- 
tion to a Stokes type differential equation of second order 
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This equation has two linearly independent solutions — the 
Airy functions of first kind A/(x) and second kind Bi{x) [15]. 
The linear combination of these two functions has to be 
normalized due to the boundary conditions. Finally we re- 
substitute Eqn. (8) to the final analytic form of l{t). For a 
more detailed description we refer to the Appendix and the 
referenced literature [16, 17]. Thus l(t) is defined as 
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It is interesting to note that the first order of the series expan- 
sion of Eqn. (10) is exactly the l{t) of the low conductivity 
solution (Eqn. (6)) as used in Ref. [11]. 

Eqn. (10) provides a complete analytical description of the 
CELIV framework. Unfortunately, the argument of the Airy 
functions in /(^) involves the parameters mobility jn and charge 
carrier density n. Therefore it is not possible to algebraically 
extract an expression for these two parameters, as it is possible 
for the mobility using a low or high conductivity approxima- 
tion [11]. However, we show later in Sec. HE that the analyti- 
cal framework can be used to derive two parametric equations 
for ju and n. Furthermore, this analytical model is capable of 
determining the experimental parameters from photo-CELIV 
experiments by fitting the photo-current transient to Eqn. (2). 
Unfortunately the Airy function of second kind Bi{x) and its 
derivative rises steeply even for small values of the argument 
X > and the fitting routine quickly gets computational can- 
cellation problems. To solve this issue we first do a scaling 
of the involved parameters to render the whole framework di- 
mensionless (see Sec. II C) and, secondly, use the confluent 
hyper-geometric function o^i [16, 18] to represent the Airy 
functions (see App. B). 
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These expressions depend on parameters which are not ex- 
perimentally accessible, such as tmax in Eqn. (17) and tmax^'^'^ 
in Eqn. (18) respectively. To overcome this issue we can use 
the scaled system and calculate the current density j{t) for a 
set of parameters. From the resulting transients we can re- 
late some of its characteristics, namely t^ax^ jo and A', 
to each other. All of these depend on A' only because of the 
scaling, which yields a definition of j{t) (Eqn. (16)) which 
is parametric in A' only. These relations between the charac- 
teristics and A' are displayed in Fig. 3 (c) (see Ref. [11, 19] 
for comparison). In Sec. HE we show how to derive para- 
metric approximations for tmax and tmax^^^'^ therefrom, which 
transform Eqn. (17) and Eqn. (18) to depend just on exper- 
imentally accessible parameters. However, we compare the 
generalized solution with the low conductivity approach of the 
CELIV framework first. 



C. Scaling to dimensionless parameters 

To yield a dimensionless system, the extraction depth l{t) 
is related to the sample thickness d and the time t is divided 
by the low conductivity case transit time ttr, Eqn. (7). This 
results in 
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Hence the scaled device thickness is J = 1 . Furthermore, the 
scaling for j and A' is set to yield a Ricatti equation (Eqn. (3)) 
and an extraction current density equation (Eqn. (2)) which is 
parametric in the dimensionless voltage slope A' only. 
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This scaling has two advantages — first it enables the com- 
putation of CELIV extraction currents over a wide range of 
parameters and second it allows us to illustrate the difference 
between the low conductivity case and the exact case in a clear 
manner (see Fig. 2). 

Note that from Eqn. (12) and Eqn. (13) the following ex- 
pression for the mobility ju and the charge carrier density n are 



D. Comparing generalized analytical and low conductivity 
CELIV framework 



In Fig. 2 we used the scaling to compare the low conduc- 
tivity and generalized case regarding the extraction depth / (?) 
and the extraction current density j{t) on a logarithmic time 
scale parametric in the dimensionless voltage slope A'. The 
black arrow indicates the direction of increasing A'. Fig. 2 (a) 
demonstrates scaled extraction current density transients j{t). 
The thin lines were calculated with the low conductivity ap- 
proach and the bold without it. 

Due to the scaling of the time axis, in case that the tran- 
sient drops at f = 1 the low conductivity approximation and 
the analytical solution are almost equal. The other transients 
show a clear distinction, with a long extraction tail for the 
general analytical solution. This long tail turns up when the 
charge carrier drift due to the external field and due to the 
redistribution of the internal field in Eqn. (2a) are balanced 
(A'/t ^ enl{tY / (2££o)). In that case the shape of the dimen- 
sionless extraction depth l{t) turns from a parabola I ocfi io 
loc VF. 

This is demonstrated in Fig. 2 (b). Here the thin dot- 
ted straight line is the extraction depth / calculated from 
Eqn. (6) (low conductivity), the buckled bold curves repre- 
sent Eqn. (11) (general solution). Hereby this straight line has 
a slope of 2 for all values of A\ The general solution qual- 
itatively follows this shape depending on A' up to a certain 
time and then deviates from the low conductivity approxima- 
tion. The quantitative deviation is plotted in Fig. 2 (c) in terms 
of the relative error of the low conductivity extraction depth 
81 (t) = 8/(^). The solid black line in Fig. 2 (b) indicates the 
scaled device thickness d. When = d all charges have 
been extracted from the device. 
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FIG. 2. Results for the calculation of (a) the dimensionless extrac- 
tion current j, (b) the dimensionless extraction depth / and (c) its 
relative error 5/ in case of the low conductivity case. The dimension- 
less voltage slope A' rises between —6.64 < log(A') < 1.36 for the 
different transients indicated by the black arrow in each of the sub- 
figures. More precise the bold dash-dotted lines corresponds to: ( - - ) 

log(AO = 1.36, ( - • - ) log(AO = -1.30, ( ) log(AO = -3.94, 

( ) log(AO = -6.64. 

(a) The thin lines are the calculated transients for the low conductiv- 
ity approximation and the bold ones are the results from the general 
scaled CELIV framework, (b) The black dotted line represents the 
low conductivity approximation and the bold ones Eqn. (10). As a 
guide to the eye the black solid horizontal line marks the dimension- 
less device thickness d. (c) The black solid line gives the relative 
extraction depth error at the time of extraction hl{ttr) and the black 
dashed line displays the relative transit time error time bttr on the 
X-axis. 



The black solid line in Fig. 2 (c) represents the relative error 
hl{ttr) at time t when l{t) = d. From the time of intersection 
tint of this black solid line with any of the bold lines the relative 
transit time error can be calculated: 5^^^ = tm — 1 • The relative 
error in transit time 5^^^ is represented by the black dotted line 
which shows 3l{3ttr) — 8/ on the y-axis versus 5^^^ on the 
X-axis. 

From this two curves we see that especially for negative 
values of log(A') the relative error 3l{ttr) gets very large, as 
well as the error in transit time 5^^^. Furthermore, we can es- 
timate the boundary of the low conductivity case validity and 
thus of the mobility equation derived from it (Eqn. (11) in 
Ref. [11]). Setting the limit for the required transit time error 
to dttr < 10% we can derive the following condition for the 
dimensionless voltage slope A' from the scaled framework 

A'>1. (19) 



For instance, in Fig. 2 (c) the second grey solid line from the 
bottom is A' = 1.12 and the time at intersection tint = 1-12. 
Hence the error is 6ttr ^ 10%. With values for common ex- 
perimental and material parameters (A' = 8 • 10^ V/s, J = 
105 nm, ^ == 6 • 10^^ m"^ and 8 = 3.7) [5, 20] A' > 1 cor- 
responds to a mobility jn < 1.69 • 10~^ cm^/Vs. 



E. Deriving parametric equations 

The dimensionless parameters in Eqn. (17) and in Eqn. (18) 
inhibit the determination of the mobility ju and the charge 
carrier density n from experiments. However, the scaled di- 
mensionless CELIV framework affords the opportunity to de- 
rive parametric equations for the charge carrier mobility ju 
and the charge carrier density n from the characteristics of 
CELIV transients, in particular from tmax and tmax^'^'^ (see 
Eqn. (17), Eqn. (18)). Whose evolution with varying dimen- 
sionless voltage slope A' is shown in Fig. 3 (c) together with 
Aj/ Jo, which equals the ratio of the real unsealed parameters 
Aj/ jo. From Fig. 3 (c) tmax and tmax^'^'^ can be related to 
Aj/ jo. Anyhow, analytical definitions for these relations are 
not available and therefore we approximate them by means 
of parametric definitions. This finally yields expressions for 
the mobility jn and the charge carrier density n including ex- 
perimental parameters only. In the following we present the 
parametric results we determined. Furthermore, we compare 
the mobility equation Eqn. (24) with parametrizations previ- 
ously suggested in literature [5, 11, 13, 19]. 

In Fig. 3 (a) the calculated relation between tmax^'^'^ and 
A j/ jo (black circles) and the parametric equation we found 
(red line) are compared. The calculated curve has a slope 
of — 1 for small values of A j/ jo < 1 and larger values 
yield a slope ^ — 2. To describe this evolution we found a 
parametrization that describes this shape with a root mean 
square deviation a = 0.8%, 

' / A -x -1.055 

tmaxA'^-^ = 0.455^ 1+0.238^ . (20) 

Aj V Jo y 

From Fig. 3 (a) it becomes clear that this approximation yields 
reasonable good fits. The relative error for all values is smaller 
than 2%. By means of Eqn. (18) and Eqn. (20) the following 
equation can be used for the determination of the charge car- 
rier density 

In case of the mobility equation Eqn. (17) parametric ap- 
proximations for tmax had been suggested in literature [5, 11, 
13, 19]. The relative errors of which are shown in Fig. 4. 
For low conductivities attended by small values of A j/ jo <C 1 , 
tmax can be analytically derived as tmax = \/l/3 [11]. How- 
ever, for common experimental conditions this equation does 
not provide adequate accuracy and it remains impossible to 
analytically define tmax in the general case. Thus, numerical 
estimated corrections have been predicted to account for the 
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FIG. 3. Calculated relation and parametric approximations for (a) 
WxV^ and (b) tmax on AJ/Jq. (c) Overview over the relation 
between the involved scaled characteristic parameters A j/ Jq (solid 
line), tfnax (dashed line), tmax^'^'^ (dash dotted line) and the scaled 
voltage slope A^ 



redistribution of the electric field. In Ref. [5, 13, 19] 
the same type 
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with a correction factor %. Juska etal [19] found % = 0.36, 
Deibel [5] suggested % = 0.21 and Bange etal. [13] published 
X = 0.18. Bange etal. suggested an additional parametriza- 
tion with a linear combination of two exponentials and four 
numerically derived adjusting parameters, which yields a very 
good fit for Aj/jo < 7 (Ref. [13] Eq. (4)). In Fig. 3 (b) 
and Fig. 4 the curves denoted by "Bange" refer to the latter 
equation from Ref. [13]. Anyhow, in Fig. 3 (b) none of the 
parametric approximations provides a good fit over the en- 
tire range of Aj/ jo. Accordingly, the best approximation we 
found giving reasonable results over the entire range of exper- 
imental parameters is 
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We clearly see in Fig. 3 (b), that Eqn. (23) renders the cal- 
culated relation very well, instead all the other parametriza- 
tion deviate sooner or later. Finally we substitute t^^^ with 
Eqn. (23) in Eqn. (17) and yield 
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for the mobility. 

In Fig. 4 we show unsealed real parameters to compare 
the different mobility equations under experimental condi- 
tions. Therefore we have calculated photo-current tran- 
sients within the CELIV framework with a defined mobility 
ju = 10~^ cm^/Vs and varied the applied extraction voltage 
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FIG. 4. Comparison of the relative mobility error 5^ for the 
different mobility equations from Juska etal. [19], Deibel [5], 
Bange etal. [13] and Eqn. (21) together with the mobilities calcu- 
lated for the low conductivity approximation from Eqn. (7) versus 
the voltage slope A^ Parameters used to calculate the transients are: 
mobility jj = 10~^ cm^/Vs, device thickness d = \50 nm, charge 



carrier density ;i = 6 • 10' 



,22 



dielectric constant 8 = 3.3. The 



dashed line is added for orientation and represents an relative error 
of 5%. The dash-dotted line is Aj/ jo. Therewith we can relate the 
relative mobility error curves to the general case in Fig. 3. The top 
axis holds the dimensionless voltage slope A' for comparison. 



slope 10^ V/s <A'< 10^ V/s. Plotted is the relative error 
5^ for the different mobility equations compared to the mo- 
bility jii we used as input parameter. To relate these values to 
the general case we added the variation of A j/ jo with varying 
voltage slope, as well as the top axes which represents the di- 
mensionless voltage slope A\ Again our parametric mobility 
equation (Eqn. (24)) is most suitable, as it yields an relative 
error which is smaller than 5% over the entire range of exper- 
imental conditions. 



III. RESULTS 

A. Fitting CELIV experiments 

In this section we briefly illustrate the computational deter- 
mination of the charge carrier transport parameters from CE- 
LIV experiments by fitting them to our analytical framework. 
Fig. 5 shows a typical photo-current transient measured with 
photo-CELIV at r = 180 K and a delay time tdeiay = 20 jus 
between laser excitation and extraction on a poly(3-hexyl 
thiophene-2,5-diyl):[6,6]-phenyl-C6i butyric acid methyl es- 
ter bulk heteroj unction solar cell. These measurements have 
already been presented by Deibel etal. [5, 20] — detailed ex- 
perimental conditions and results can be found there. The pa- 
rameters in Eqn. (2) to be fitted to the experimental signal 
are the mobility ju and the charge carrier density n. The fit- 
ting results are summarized in Tab. I together with the values 
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FIG. 5. Results fitting photo-current transient from a photo-CELIV 
measurement aiT = 180K and a delay time t^ieiay = 20jus. The black 
circles connected with a black thin dashed line correspond to the 
measurement. The thick dashed line represents the fit with a unipo- 
lar extraction current. For the solid line an ambipolar extraction is 
assumed for the fit. The two dash-dotted lines are the hole and the 
electron extracting currents, which superpose to the resulting photo- 
current (solid line). 



determined from the measurement with Eqn. (24) (mobility 
ju) and Eqn. (21) (charge carrier density n). We find that fit- 
ting with one extraction current (unipolar transport), taking 
one conducting type of charge carrier into account, fails, see 
Fig. 5 (long dashed line). Here the fitted charge carrier den- 
sity is too high and the time of the extraction current max- 
imum tfnax is too short, thus the charge carrier mobility is 
overestimated. If we instead assume an ambipolar extraction 
to fit the photo-current, the resulting curve in Fig. 5 (solid 
line) matches very well with the measured data. For the am- 
bipolar extraction we used, as a first order approximation, a 
linear superposition of two extraction currents. The nature 
of the latter cannot be determined experimentally, but due to 
the photo-generation of excitons by laser excitation it is very 
likely to be hole and electron driven currents, respectively. 
As we see from the dash dotted and the dash double dotted 
lines in Fig. 5 the extraction current maxima are only slightly 
shifted with respect to each other, implying that the mobil- 
ities of holes jLihjit and electrons jHeju are almost balanced, 
see Tab. 1. This is in a good agreement with time-of-flight 
measurements published in Ref. [7]. Correspondingly, we as- 
sign the slightly slower charge type as hole and the other one 
as electron transport to simplify the discussion. Furthermore 
the electron densities rieju are about a factor of two smaller 
compared to the hole densities Uhju implying that the electron 
density is more strongly reduced within the delay time t delay 
Therefore we propose the following two explanations with- 
out going too much into detail. Firstly, enhanced trapping of 
the electrons before the extraction seems possible. However, 
PCBM is found to be a trap-free acceptor [21-23], thus the 
trapping could not be energetically, but more a morphological 
effect, where the electrons are immobilized in isolated phases. 
For the small spherical fullerenes [24] this is more likely than 



for the long polymer-chains. Secondly, in the photo-CELIV 
method an offset voltage is applied which compensates the 
built-in field to prohibit the extraction of charges before be- 
ing swept out. However, the recombination of charge carriers 
changes the flat-band conditions continuously with time, thus 
a certain amount of charges is extracted during the delay time. 
Hence, electrons in conducting phases are more efficiently ex- 
tracted due to their slightly higher mobility than holes. 

Finally, we combine in Tab. I our fitting results (jnhjit, l^ejit^ 
l^unu ^hjit^ ^ejit^ ^uni) and the valucs we derived from the mea- 
surements with the help of Eqn. (24) (jHcaic) and Eqn. (21) 
(ricaic)- In addition, we define a weighted mean mobility 
Jiimean = {nhjitl^hjit ^nejiti^ejit) / {nujit ^Uejit) as wcU as the 
total charge carrier density Utotal = ^hjit + ^ejit from the am- 
bipolar fit to compare it with the calculated values jUcaic and 
ritotal- Furthermore, the charge carrier density is determined 
from the area below the photo-current transient (riint), as it is 
commonly done [5, 20, 25, 26]. 

From the photo-CELIV measurement we obtain a fitted 
hole mobility jUhju comparable to the calculated jHcaic- In- 
stead, the weighted mean mobility jUmean is slightly above 
jLicaic and the fitted electron mobility jHejit is two times higher 
than jLicaic- Thus, we assign the mobility determined from the 
photo-current peak maximum to the mobility of the holes, due 
to the higher hole density (see Tab. I). In general the mo- 
bility JLicaic is related to the more conducting charge carrier 
type [27] . From our results we conclude, that the the charge 
carrier density is the crucial parameter determining the charge 
carrier type, which is related to the mobility jUcaic- Moreover, 
the extraction current (Eqn. (2)) depends quadratically on the 
charge carrier density and only linearly on the mobility. 

We note that the integrated charge carrier density Hint ac- 
counts for the density of extracted charges, while the values 
from the transient fitting and from Eqn. (21) reflect the total 
density of mobile charges involved in the photo-current. As 
expected the density of extracted charges Hint is smaller than 
the values derived from Eqn. (21) ricaic and from the ambipo- 
lar fit ritotal • About 65 % of the photo-generated charge car- 
riers could be extracted within the length of the applied volt- 
age pulse tp = I ms. Experimentally, this limitation is often 
observed, when the trade-off between the applied maximum 
voltage, the signal to noise ratio and the extraction time pre- 
vents the complete extraction of all free charge carriers. The 
tail slope of the extraction current and hence the amount of ex- 
tracted charges gets very small. Thus even a two times longer 
extraction time is not capable to extract all charges. However, 
the total density of mobile charges derived from the ambipo- 
lar fitting ritotal and from the calculation via Eqn. (21) ricaic are 
in a good agreement. Therefore and due to small root mean 
square deviation in Fig. 3 (a), we recommend using Eqn. (21) 
to determine the charge carrier density n from experiments. 

We want to point out that ricaic should to be the upper limit 
of extracted charges. If this is not the case, the photo-current 
transients are distorted by processes which are not consid- 
ered in the general CELIV framework, for instance recombi- 
nation [13] or trapping. In such a case the evaluation with any 
equation derived there from could be questioned. However, as 
shown in our measurements the extracted charge carrier den- 
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ambipolar 


uiiipoidi 






holes electrons mean/total 




integrated this work 








^uni 


^int ^calc 


mobility ju 
(xlO--^^) 


1.52 


3.09 2.30 


3.69 


— 1.64 


density n 
(xlO-l^cm-^) 


2.67 


1.08 3.75 


7.68 


1.67 4.65 



TABLE I. Summary of the extracted parameters for the photo-CELlV measurement shown in Fig. 5. For the sake of clarity, the mobihty /i 
and the charge carrier density n is represented by the x in the table header. The parameters are taken from the unipolar {xuni) and ambipolar 
(xhjit,Xejit) fits or calculated {Xcaic)- The mobility is calculated from Eqn. (24), the charge carrier density from Eqn. (21) and by integrating the 
extraction current density {nint)- The ambipolar weighted mean mobility is approximated by fimean = (p^hjitl^hjit ^^ejitl^ejit) / (p^hjit ^ejit) 
and the ambipolar total charge carrier density is Utotal = ^hjit ~^^ejit- 



sity Hint is smaller than the calculated Ucaic 



IV. CONCLUSION 

With the analytical solution for the extraction depth l{t) 
we derived a complete framework for the CELIV technique. 
Therewith the extraction current response due to a linearly in- 
creasing voltage can be analytically calculated. 

We suggested two new parametric equations for the deter- 
mination of the charge carrier mobility ^ and the charge car- 
rier density n from the characteristics of CELIV transients, 
respectively. These equations are capable of handling the en- 
tire experimental range of parameters. The relative error of 
the charge carrier density equation does not exceed 2% and 
accounts for the total density of mobile charges and not just 
for the extracted charge carrier density. Our mobility equation 
yields lower deviations from the analytical predictions than 
any previous suggested mobility equation. 

Finally we evaluated photo-CELIV measurements by fit- 
ting them within the derived analytical CELIV framework. 
We find that reasonably good fits can only be achieved, if an 
ambipolar extraction of holes and electrons is taken into ac- 
count. The results show a balanced hole and electron mobility 
in P3HT:PCBM solar cells in accordance with previous ex- 
periments [7]. Furthermore we found that the type of charges 
with the higher charge carrier density, instead of the more mo- 
bile one, is mainly rendering the shape of the photo-CELIV 
transients and is therefore the one that is probed. 
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Appendix A: The Airy functions 

We briefly want to show how to normalize the Airy func- 
tions to yield the solution of x(t) Eqn. (9). The Airy functions 
are defined for real values x as follows 




and the values of Ai{x) and Bi{x) at x = are 

^'(°) = 3^(273) ^^^^ 

The solution for Eqn. (9) is the linear combination of Ai{x) 
and Bi{x) 

x{t)=CiAi'{x)^C2Bi'{x) , (A5) 

We can derive the following boundary conditions forx(O) and 
jc'(O) from /(O) = at time t = 

jc(0) = const. (A7) 
y(0)=0. (A8) 

Therefrom we can determine Ci and C2 to 

Ci = ^3^/^r(2/3) (A9) 
C2 = ^3i/6r(2/3) . (AlO) 

Appendix B: Alternative definition of the extraction depth 

The Airy functions in Eqn. (10) can be revealed by the con- 
fluent hyper-geometric function 0^1 • This yields a computa- 
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tional more robust definition of the extraction depth. 
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